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Experimental vibrational spectra (Raman and infrared absorption) 
of berberine are obtained at room temperature. The vibrational 
spectra of berberine are calculated by the DFT method at the 
B3LYP/6-31H — hG(d,p) level. Based on the correlation between 
experimental and calculated data, the vibrational spectrum is in- 
terpreted in the frequency range of 800-1700 cm -1 in detail. The 
experimental and calculated spectra of intramolecular vibrations 
are found to correlate closely. 



1. Introduction 

Natural alkaloid berberine has been used in medicine for 
a long time for treating a number of diseases. Prepa- 
rations of berberine have been found to have cytotoxic, 
bactericidal, and antiviral effects [1]. Moreover, berber- 
ine influences the cancer cell metabolism by destroying 
cells irreversibly. Such activity of the alkaloid is at- 
tributed to its ability to intercalate into the DNA macro- 
molecule by blocking the processes of replication and 
transcription. Many articles are devoted to the interac- 
tion between berberine and nuclei acids, see, e.g., reviews 
[1,2]. Another important target for antitumour berber- 
ine agents is DNA-topoisomerase: berberine can inhibit 
topoisomerase, by breaking the connection between this 
enzyme and DNA [3]. 

To clarify which mechanism causes berberine to pro- 
duce its therapeutic effect, it is necessary to know the 
mechanisms by which it binds to DNA. Earlier, using the 
Raman spectroscopy method, we revealed an interesting 
fact of a resonance interaction of berberine and DNA vi- 
brations [4, 5]. It was found that, in the DNA-berberine 
solution, the intensity of vibrations of both berberine 
and DNA increases greatly (by orders as compared to 
those of the spectra of components) in the range where 



their intense spectra overlap (1000-1700 cm^ 1 ). To ex- 
plain this interaction, the detailed data about berberine 
and DNA vibrational spectra are required. 

The vibrational spectra of berberine were previously 
investigated with the use of the Raman [6-9] and IR- 
absorption spectroscopy methods [7, 9], but no exhaus- 
tive interpretation of the obtained spectra was given 
in those works. In works [7-9], the Raman spec- 
tra of berberine in the range of 600-1800 cm -1 ob- 
tained by the SERS and SSRS (surface-enhanced and 
shifted-subtracted Raman spectroscopy) methods were 
presented. 

This work continues and develop our study described 
in [10]. 

2. Materials and Methods 

Natural berberine is mainly contained in plants as 
hydrochloride or dihydrosulfate, a structural formula 
of a neutral berberine molecule is C20H19NO5 (or 
C2oHi8N04(OH)). A structural formula of a berberine 
cation (as usually investigated) is [C2oHi§N04] + (see 
Fig. 1). The quantum- mechanical investigations of the 
berberine structure were carried out in [11, 12]. The 
results of studies of the spatial structure of berber- 
ine cation isomers [12] demonstrate that the berberine 
cation that consists of hexamerous planar rings has an 
almost planar frame structure. A deviation from this 
planar structure is observed only in a partially saturated 
C ring. It conforms to the shape of a half-chair, in which 
C7 and C8 atoms of the C ring deviate appreciably from 
the plane of B and D rings. 

The microcrystalline powder of berberine has been 
studied (berberine hydrochloride, "Alps Pharmaceuti- 
cal", Japan, 407 g/Mole). Since berberine absorbs in the 
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range to 550 nm [13], a He-Ne laser emission (A = 6328 
A) has been used for the excitation of Raman spectra; 
the real power on a specimen was ~ 10 mW. Some spec- 
tra were obtained by using a Kr-laser (A = 6471 A) with 
a power ~ 60-80 mW on a specimen. Raman spectra 
have been measured at room temperature with the 45°- 
reflection geometry using an optical setup made on the 
basis of a double-grating monochromator DFS-24 with 
a resolution of 1.5 em . A laser beam was focused on a 
specimen with a cylindrical lens to obtain an elongated 
irradiated area parallel to the entrance slit of the spec- 
trometer. To improve the signal-to-noise ratio, the time 
accumulation at the point was ~ 1.6 s (multipass mode). 
Spectra were recorded in the range of 40-4000 cm -1 . 

The IR-absorption spectra of berberine (microcrys- 
talline powder) were recorded with a Nicolet NEXUS- 
470 Fourier spectrometer. We used the ATR (Attenuated 
Total Reflectance) technique. To improve the signal-to- 
noise ratio, a signal accumulation regime was used (128 
scans). Spectra were recorded in the range of 400-4000 
cm with a resolution < 1 cm . In the resulting spec- 
tra, Fresnel losses on the input and output surfaces of 
specimens were allowed for. In addition, to obtain the 
absorption spectra, the dependence of the penetration 
depth of evanescent waves into the sample on their wave- 
length is considered according to the relation [14] 

d P = 1 X =, (!) 

2Trn\i sin 2 8 — (n/n s ) 2 

where n and n s are refractive indices of a prism (dia- 
mond) and a specimen, respectively, A is a light wave- 
length, and 9 is an incident beam angle. 

In some cases, experimental spectra were treated, by 
using the Origin and PeakFit programs to specify the 
data on the band position, the number of band compo- 
nents, and so on. 

The optimized geometry of the berberine cation 
was calculated by the DFT method at the B3LYP/6- 
311++G(d,p) level without any structural restrictions. 
At this level, the vibrational (Raman and IR-absorption) 
spectra were calculated in the harmonic approximation. 
It should be noted that, for the first time in an optimized 
structure, the intermolecular hydrogen bond CH. . . O 
between metoxylic groups OCH3 of berberine was found 
[10]. The existence of the bond was determined based 
on the presence of a critical point of (3,-1) type on the 
electron density distribution by the AIM method [15]. 
The binding energy was found to be 3.12 kcal/mole. 
The DFT calculations were performed using the " Gaus- 




O N 

Fig. 1. Structural formula of a berberine cation (calculation by 
the DFT method at the B3LYP/6-311++G(d,p) level) 

sianOS" program package for Win32 [16] granted by the 
" Gaussian" corporation. 

Because the DFT method usually overestimates fre- 
quencies, a scaling factor of 0.985 was used for the com- 
parison of calculated and experimental data. This cor- 
rection was necessary due to errors in the calculation 
of interatomic interactions as a result of the limited 
basic functions set [17]. The value 0.985 provides the 
best correspondence between calculated and experimen- 
tal data. The scaling factor was determined by Ra- 
man spectra - as it turned out, the shapes of calcu- 
lated and experimental spectra was more similar for Ra- 
man spectra, and the correspondence between calculated 
and experimental frequencies was practically unambigu- 
ous. In our case, the scaling factor is close to 1 indicat- 
ing that the used basic set is sufficient for the studied 
molecule. 

In addition, it is known that the Gaussian set of pro- 
grams deals with Raman activities Si, not intensities Ii 
[18, 19]. Therefore for comparing with the experiment, 
the Raman activities Si were corrected by the following 
relationship to obtain Raman intensities [18]: 

j = C(v - v. i ) i S l ^ 
4 Ui[l — exp(— hcfi/kT)] 

Here, C is a constant, vq is the laser excitation line 
frequency, Vi is the vibrational frequency, h, k, c, T are 
Planck and Boltzmann constants, speed of light, and 
temperature in degrees Kelvin, respectively. 

The values of frequencies and relative intensities of 
spectral lines with regard for the aforementioned correc- 
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2. Raman spectra of berberine in the range of 1000—1700 
1 : experimental, A cx = 6328 A (a), and calculated (6) 



tions and the detailed interpretation of vibrations are 
presented in Table. 



3. Results and Discussion 

The berberine crystal (C20H18NO4 C1~4H 2 0) symmetry 
is triclinic (PI, z — 2) [20]. Due to a low symmetry, there 
are no degenerated vibrations in vibrational spectra of 
berberine. The correspondence between the "crystalline" 
and "molecular" modes is unambiguous. The most essen- 
tial difference between the spectra of free and crystalline 
berberine had to reveal itself in the range of very low 
frequencies (the so-called external vibrations, there are 
three of them at z = 2. In the range of middle and 
high frequencies (internal vibrations), the influence of a 
crystal structure is very weak. Both duplication in the 
number of vibrations (because of z = 2) and some shift- 
ing of the vibrational bands take place. 
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Fig. 3. IR- absorption spectrum of berberine in the range of 800— 
1700 cm -1 : experimental (a), and calculated (6) 



Therefore, we believe that the comparison of experi- 
mental vibrational spectra of microcrystalline berberine 
with calculated spectra of the berberine cation in the 
actual range of 600-1800 cm -1 was sufficiently correct, 
which was confirmed by the following analysis. Experi- 
mental and calculated Raman and IR-absorption spectra 
in the range of 800-1700 cm -1 are presented in Figs. 2 
and 3, respectively. The frequency values - both ex- 
perimental and calculated using the scaling factor - are 
presented in Table. The relative intensities of Raman 
lines were corrected by (1). In this case, we are inter- 
ested more in the range of 1000-1700 cm -1 . This is be- 
cause berberine vibrations are the most intense in this 
range, and there a resonance interaction of the berberine 
and DNA vibrations (increasing by orders) is observed in 
the Raman spectrum of a berberine-DNA water solution 
[4, 5]. 
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Experimental and calculated frequency values, as well as the interpretation of vibrations. The frequency scaling factor 
is 0.985, calculated Raman intensities were corrected by relationship (2), experimental IR-intensities — by (1) 



Experimental Data 


Calculations 
with corrections 


Interpretation of vibration modes 
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strong i/(N-CH 2 ), 5(C5NC18); 
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Experimental and calculated frequency values, as well as the interpretation of vibrations. The frequency scaling 
factor is 0.985, calculated Raman intensities were corrected by relationship (2), experimental IR-intensities - by (1) 
(Continuation) 



Experimental Data 


Calculations 
with corrections 


Interpretation of vibration modes 
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Experimental and calculated frequency values as well as interpretation of vibrations. The frequency scaling fac- 
tor is 0.985, calculated Raman intensities were corrected by relationship (2), experimental IR-intensities - by (1) 
(Continuation) 



Experimental Data 


Calculations 
with corrections 


Interpretation of vibration modes 
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1609.5 


(0.4) 


(4.5) 


in-plane u-S vibrations of atoms C,0,H,N of all rings; 


sh 
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1618 


(1.5) 


(0.3) 


weak u,p(CH 3 ); r(CH 2 ) r.C, f(CH 2 ) r.E; 
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1622 


1629 


1624 


1620 


1625 


(10.7) 


(3.3) 


1609: strong of all rings; O - almost immobile; 


















1618: strong vibrations r. A,D; N - immobile; 


















1625: strong vibrations - r.D, N - vibrates; 


sh 
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1633 


1634 


(0.9) 


(0.5) 


in-plane v — <5 of atoms C,N,0,H of all rings 


















(strong - r.B,C,D, weak - gr. CH2; gr. OCH3 - immobile); 



Notes: In brackets near frequency values are the relative intensity of lines. (The intensities of the modes at 1397 cm -1 in Raman 
spectra and 1367 cur 1 in IR-absorption spectra count as 10. Modes with relative intensity less than 0.1 are designated as "vw", for 
"very weak"). 

Abbreviations: b: band; sh: band shoulder; r.A,B,C,D,E: a, B, C, D and E rings, respectively; weak" and "strong" denote 
vibrations with small and large displacements, respectively. Vibrations: v - stretching, <5 - deformation (for the CHVgroups, the same 
as bending); r: - torsion; v — <5: a part of C atoms in a ring performs i/-vibrations, and the others perform 5-vibrations; for 
CHVgroups - /: wagging; p: rocking (atoms vibrate in one plane); for CHVgroups - u: umbrella (symmetric <5- vibrations) ; p: rocking 
(atoms vibrate in parallel planes); for r, /, p - vibrations preserve the rigidity of CH3 and CH2 groups. 



Unfortunately, the high frequency region of 3000-3300 
cm -1 was not applicable to the comparative analysis: 
Raman spectra have very small intensities in this re- 
gion, and IR-absorption spectra are only loosely corre- 
lated with calculations. 

The calculated spectrum consists of 123(3 x 43 — 
6) non-degenerate vibrations of a berberine cation 
[C 2 oHi 8 N0 4 ] + , 105 of which lie in the range of 20-1700 
cm -1 (61 of which are in the actual range of 800-1700 
cm -1 ), other 18 ones fall in the range of 3000-3300 cm -1 . 
Low-range frequencies up to 720 cm -1 correspond to 
out-of-plane vibrations of rings and associated groups. 
In-plane vibrations of rings appear beginning with 720 
cm -1 ; in particular, the only sufficiently intense 727- 
cm" 1 Raman mode outside of the range of 1200-1700 
cm -1 . In-plane vibrations are characteristic of the entire 
range from 720 to 1700 cm -1 . The high-frequency region 
of 3000-3300 cm" 1 consists of the modes corresponding 
to vibrations of C-H-bonds with weak displacements of 
other atoms, mostly C. 

In Fig. 2, a, the experimental Raman spectrum of 
berberine (microcrystalline powder) in the range of 
1000-1700 cm" 1 at the excitation by the 6328-A line 
is presented. We adjusted the spectra, by subtracting 
the linearly decreased background and by smoothing 
the noise level slightly. In the experimental spectrum, 
more than twenty vibration modes were registered, and 
the strong correlation was observed between our data 
and data of [6, 9] (see Table), where the reasonably 
good berberine Raman spectra in the range of 1000- 



1700 cm" 1 were presented (in [9], FT- Raman). Unfor- 
tunately, in Tables in [6, 9], only the frequencies of the 
most intense bands were given; an interpretation of vi- 
brations was absent, or it was incomplete or inadequate. 
It should be noted that the Raman spectra of berberine 
obtained in [7-9] by the SERS and SSRS methods are 
not applicable fully for the correct comparative analysis 
of the experimental and calculated data. 

In Fig. 2,b, we present the calculated Raman spectrum 
of a berberine cation in the range of 1000-1700 cm" 1 , 
by using a scaling factor. 

The frequency values and the detailed interpretation 
are presented in Table. The strong correlation between 
experimental and calculated spectra is observed in the 
range of 1000-1700 cm" 1 - both between frequencies and 
intensities of vibrations. In the range up to 1000 cm" 1 , 
besides a 727 cm" 1 mode, the intensities of calculated 
vibrations are very small (less by 2-4 orders) in agree- 
ment with experimental data. Notice that some modes 
very weak by calculations (e.g., 1302 cm" 1 ) are revealed 
in experimental spectra. 

In Fig. 3, a, the IR-absorption spectrum of berber- 
ine (microcrystalline powder) in the range of 800-1700 
cm" 1 is presented (after subtracting the background), 
about 40 rather intense separated vibrational modes are 
observed, the good agreement between our data and data 
of [7-9] takes place (see Table). 

It should be mentioned that the obtained experimental 
IR-absorption and Raman spectra appear quite similar 
to the IR-absorption spectra and the FT-Raman spec- 
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trum presented in [7, 9]; but IR- frequency values were 
presented in [7] only for the most intense lines, whereas 
they were not presented at all in [9]. Moreover, the anal- 
ysis in [7, 9] was performed with SERS spectra that did 
not correspond well to our spectra by shape. The cal- 
culated Raman and IR-absorption spectra appear very 
similar to those we calculated, but differences in fre- 
quency values are quite substantial. Experimental and 
calculated Raman spectra are similar, but IR spectra are 
not. 

We attribute this discrepancy to the fact that the cal- 
culations usually use an ion or molecule of berberine, 
but experiments are conducted with commercially pro- 
duced crystallohydrate. Since the dipolarity of vibra- 
tions for IR-absorption processes is important, we have 
not discounted the possibility that the presence of dipole 
water molecules (and OH-groups) in specimens might 
have led to the dipole-dipole interaction with berberine 
molecules, and the results of the interaction might be 
revealed in IR-absorption spectra. The influence on Ra- 
man spectra is weaker, because the Raman scattering 
processes are related to the electron system, and the in- 
fluence of the dipole-dipole interaction on the vibrational 
system is less direct. 

In Fig. 3,6, we present the calculated IR-absorption 
spectrum of berberine (cation) in the range of 800-1700 
cm , by considering a scaling factor. 

The relatively strong correlation between the frequen- 
cies of experimental and calculated IR-absorption spec- 
tra is observed in the range of 800-1700 cm -1 ; the cor- 
relation between the intensities of vibrations is some- 
what worse, for the aforementioned reasons. It should 
be noted that, contrary to Raman spectra, some modes 
that were calculated to be intense were not observed at 
all in experimental IR-spectra, or they had only a very 
weak intensity (e.g., 1149, 1242, 1285, 1408, 1471 cm" 1 , 
etc.). In this case, the lines of an unknown nature were 
revealed in Raman and IR spectra (1127 cm -1 (IR), 1129 
cm -1 (Raman) etc.). 

4. Conclusions 

Thus, the Raman and IR-absorption spectra have been 
obtained for a berberine cation by DFT at the B3LYP/6- 
311++G(d,p) level. The optimized geometry of the 
berberine cation was calculated as well. The results of 
calculations strongly correlate with experimental data 
obtained for microcrystalline berberine chloride in the 
region of 1000-1700 cm -1 significant from the viewpoint 
of the interaction of berberine with DNA. The obtained 
interpretation of the Raman bands of berberine can be 



used for an analysis of the berberine interaction with 
nuclei acids and other biomolecules. The DFT method 
can be used, with satisfactory reliability, to calculate vi- 
brational spectra of other isoquinoline alkaloids of the 
protoberberine group, for which obtaining the Raman 
vibrational spectra would be complicated or even im- 
possible for one reason or another. 

The autors are grateful to S.O. Alekseev and 
M.E. Kornienko for the help in the recording of IR- 
spectra and to L.A. Zaika (Institute of Molecular Biology 
and Genetics of the NANU) for specimens of berberine. 
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KOJIHBAJIBHHH CIIEKTP OPrAHIHHOI CnOJIYKH 
BEPBEPHHY TA HOrO IHTEPIIPETAIILH 
KBAHTOBO-MEXAHIHHHM METQUOM 
OyHKHIOHAJIA ryCTHHH 

H.B. BauiMCiKoea, CIO. Kymoeuu, P.O. XypaKiecbKuu, 
fl.M. roeopyn, B.M. H-w,yK 

P e 3 io m e 

3a KiMHaTHoi' TeivinepaTypn OTpuMaHO KOjinBajibm cneKTpn (pa- 

M&HlBCbKIIH Ta iHfppaMepBOHOrO nOrjIHHaHHJl) MlKpOKpiICTajli^IHO- 

ro xjiopHfly 6ep6epHHy Ta npoBefleHO ixhio iHTepnpeTaii,iio Me- 
TOflOM <py HKD ,i° Hajia ryCTHHH Ha piBHi Teopi'i DFT B3LYP/6- 
311++G(d,p) b fliana3om nacTOT 800-1700 cm -1 . CnocTepe»ceHO 
flo6py KopejiHr;iio Mi>K eKcnepiiMeHTajibHHMH Ta po3paxoBaHHMn 

x iaCTOTaMH KOJIHBaHB. 
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